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ABSTRACT 


The  instantaneous  surface  regression  rates  have 
been  calculated  as  functions  of  time  and  axial 
distance  from  the  radius  vs  time  data  obtained 
from  plexiglas  combustion.  Temperature  measure¬ 
ments  of  the  fuel  behind  the  regressing  surface 
have  been  made.  Additional  experiments  were 
carried  out  to  examine  the  burning  rates  of  dif¬ 
ferent  fuels  such  as  polystyrene,  polyethylene, 
30  and  60  ^  aluminum,  30  ^  titanium  hydride,  and 
1/4  and  1/2  ^  carbon  in  plexiglass, ^^S^e  sig¬ 
nificance  of  the  burning  rate  and  temperature 
measurement,  results  is  discussed. 
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INTRODUCTION  AND  OBJECTIVES 


This  research  program  in  hybrid  combustion  has  been  concerned  with  the 
investigation  of  the  processes  (both  physical  and  chemical)  involved  in 
the  interactions  between  solid  polymeric  fuels  and  gaseous  oxygen.  The 
objective  of  these  studies  has  been  to  determine  experimentally 
which  processes  are  important  in  affecting  the  surface  regression  rates i 
and  to  develop  a  suitable  theoretical  model  to  account  for  the  observed 
results . 

A  simplified  theory  of  hybrid  combustion  as  sophisticated  as  the  classi¬ 
cal  treatment  of  Burke  and  Schumann  (Ref.  1  )  for  the  flame  shape  in 
cylindrical  diffusion  flames  was  presented  for  the  plexiglass-oxygen 
(  system  (Ref.  2  ).  This  theory  involved  a  model  obtained  using  the 

following  assumptions: 

1.  Chemical  reactions  are  confined  to  the  flame  front,  considered 
an  infinitely  thin  mathematical  surface  close  to  the  wall,  to 
which  the  fuel  and  oxidizer  are  delivered  in  stoichiometric 
proportions  and  are,  therefore,  consumed  instantaneously.  For 
this  reason,  no  fuel  exists  in  the  core  and  no  oxygen  exists 
in  the  annular  region  between  the  flame  front  and  the  wall. 

2.  The  pressure  gradient  is  negligible  and  the  linear  flow  veloci¬ 
ties  are  constant.  Hence,  the  axial  velocity  term  is  neglected, 
and  only  radial  transport  by  diffusion  is  considered.  This  has 
been  assumed  to  be,  primarily,  concentration  diffusion. 
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3.  Gas  mixtures  are  considered  to  be  binary  mixtures  composed  of 
the  major  constituent  (either  fuel  or  oxidizer,  depending  on 
location  beneath  or  above  the/ £^e  frx>ni)  and  the  rest  of  the  gas 
mixture  as  an  average. 

4.  The  product  of  gas  density  and  binary  diffusion  coefficient  is 
assumed  to  be  independent  of  radial  and  axial  distances  and 
temperature. 

3.  No  distinction  is  made  between  types  of  flow,  and  appropriate 
applicable  transport  coefficients  are  to  be  used  for  the  axi- 
symme trie  c ond i t i ons . 

6.  The  binary  diffusion  coefficients  of  fuel  and  oxygen  through 
nonreacting  gas  are  numerically  equal. 

With  proper  definition  of  the  diffusion  parameters,  this  theory  may  be 
applied  to  all  types  of  flow  and,  in  particular,  to  the  turbulent  flow 
conditions  examined  experimentally. 

The  purpose  of  the  more  recent  work  has  been  to  define  better  the  contri¬ 
butions  of  the  experimental  parameters  to  the  observed  regression  rates 
(Ref.  3  through  5  )  and  to  determine  the  validity  of  the  above  assump¬ 
tions.  In  this  reporting  period,  effort  was  directed  toward  obtaining 
better  defined  functions  of  radius  with  time  at  the  two  entrance  Reynolds 
numbers  (55,000  and  36,000)  and,  from  these  functions,  calculating 
instantaneous  regression  rates.  In  addition,  the  temperature  profile 
in  the  solid  fuel  ahead  of  the  regressing  surface  was  further  examined, 
and  the  isothermal  polymer  pyrolysis  rates  and  surface  regression  rates 
were  compared . 
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EXPERIMENTAL 


The  basic  apparatus  has  been  described  in  previous  reports  (Ref.  2  and 
3).  The  tubular  combustor  configuration  urats  used  for  the  experimental 
program.  The  mixing  and  flov  stabilization  tube,  vhich  vas  a  2-inch 
square  duct,  has  been  replaced  by  a  1-inch  diameter  tube  to  reduce  the 
turbulence  caused  by  the  abrupt  change  in  diameter  at  the  entrance.  In 
a  typical  run,  a  fuel  cylindrical  surface  composed  of  3  by  3  by  2-inch 
blocks  with  1-inch  diameter  concentric  holes  aligned  with  the  holes  of  the 
facing  and  backing  plates,  is  burned  for  the  specified  time.  The  weight 
lost  by  each  block  is  measured  to  determine  the  surface  regression;  axi- 
symmetric  burning  is  assumed.  The  axial  length  (niunber  of  blocks),  the 
composition  and  temperature  of  fuel,  the  time  of  burning,  and  the  initial 
diameter  of  the  concentric  hole  may  be  varied.  Thermocouples  may  be  im¬ 
bedded  between  adjacent  fuel  blocks  at  any  desired  radial  distance  to 
measure  the  increase  in  temperature  at  that  position  as  the  burning  sur¬ 
face  approaches.  The  latest  thermocouples  used  in  this  study  were  chromel- 
alumel  of  either  1-  or  0.6-mil  diameter  wire. 

The  thermogravimetric  balance  apparatus  has  been  described  in  Ref.  2. 

The  experimental  procedure  consisted  of  observing  the  weight  loss  of  a 
block  of  polymer  as  a  function  of  time,  at  a  constant  temperature,  and 
in  the  presence  of  different  oxygen-nitrogen  mixtures. 
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RESULTS 

TIME  DETENDEMCE  OF  RADIUS  AND  BURNING  RATES 

To  understand  the  hybrid  burning  process  better,  it  vas  believed  that  in¬ 
stantaneous  burning  rates  at  definite  axial  positions  should  be  examined, 
since  these  would  be  more  meaningful  than  time  and  length  dependent 
averages.  To  do  this,  it  was  necessary  to  make  runs  at  several  values  of 
t  and  then  derive  an  r(t)  equation.  This  was  done  with  plexiglass 
as  a  fuel,  primarily  at  two  entrance  Reynolds  numbers,  35>000  (22.0  gm 
02/sec)  and  36,000  (l4.4  gm  02/sec),  with  additional  experiments  at  6OOO 
(2.4  gm  02/360).  On  ignition,  a  line  pressure  of  about  3  or  4  psig  was 
observed  at  the  55 1 000  flow  which  reduced  to  zero  after  about  10  to  15 
seconds.  The  lower  flowrates  exhibited  no  pressure  rise  in  the  line. 

It  is  believed  that  the  small  pressure  rise  would  have  a  negligible 
effect  on  the  burning  rate. 

By  use  of  the  Research  Multiple  Correlation  Computer  Program,  \diich  applies 
various  statistical  tests  attempting  to  obtain  the  best  fit  of  the  supplied 
functions  to  the  experimental  data,  Eq.  1  was  obtained  for  radius  vs 
time: 


r  =  A  -I-  Bt 


0.8 


+  Ct 


0.3 


(1) 


where 


A,  B,  and  C  -  constants 
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Table  1  summarizes  the  data,  and  Fig,  1  and  2,  which  present  a  plot 
of  final  radius  vs  time,  illustrate  typical  fits  of  data  and  equation. 
The  same  equation  was  used  for  all  blocks  at  the  higher  flowrates  and 
for  the  entire  time  range.  It  must  be  re-emphasized  that  £q.  1  was  not 
intended  to  have  theoretical  significance.  It  was  to  be  used  only  to 
determine  the  time  and  axial  dependences  of  the  regression  rates. 

The  curves  presented  in  Fig.  1  and  2  are  not  very  accurate  at  times 
less  than  10  seconds,  and  the  slopes  go  to  infinity  at  t  -  0.  A  perfect 
fit  would  require  a  sigmoidal  curve  to  indicate  the  surface  heatup  time 
just  following  ignition.  According  to  the  data,  t  for  inflection  of  the 
sigmoid  (point  of  maximum  slope)  would  be  below  10  seconds,  but  accurate 
data  could  not  be  obtained  for  these  short  times.  Thus,  the  heatup  time 
was  not  defined  experimentally. 

By  differentiation  of  Eq.  1,  the  regression  rate  as  a  function  of  time 
can  be  obtained: 

r  =  0.8  +  0.3  (2) 

The  results  of  these  calculations  are  summarized  in  Table  2.  Figure 
3  and  4  illustrate  typical  examples.  It  is  interesting  to  note  that 
the  regression  rates  at  the  6000  flow  are  about  constant  with  time  and 
axial  distance,  although  the  average  rates  (over  time  and  distance)  de¬ 
crease  somewhat  with  time,  as  indicated  on  the  following  page. 
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TAB1£  1 

CONSTAOTS  FOR  RADIUS  VS  TIME  EQUATION 
r  =  A  -h  (cm) 


Oxygen 

Flov 

36,000,  21  Runs 

55 9 000,  22  Runs 

BlociV^on  stan  1 8 

A 

B 

C 

0^ 

A 

B 

c 

(T» 

XlOO 

XIO 

XlOO 

XlOO 

XIO 

XlOO 

1 

0.877 

0.183 

2.59 

1.5 

0.940 

0.176 

2.42 

1.2 

2 

1.040 

1.16 

1.25 

1.2 

1.101 

1.38 

1.03 

1.2 

3 

1.131 

1.31 

00 

0 

1.0 

1.186 

1.61 

0.550 

1.0 

4 

1.139 

1.24 

0.782 

1.2 

1.160 

1.42 

0.786 

0.9 

5 

1.086 

1.12 

1.11 

1.5 

1.083 

1.18 

1.27 

1.0 

6 

0.993 

0.98 

1.53 

1.1 

0.992 

0.933 

1.80 

1.1 

7 

0.976 

0.98 

1.58 

1.5 

0.987 

0.947 

1.84 

1.5 

8 

1.000 

1.10 

1.41 

1.5 

1.021 

1.07 

1.70 

2.4 

9 

1.090 

1.42 

0.896 

1.5 

1.029 

1.16 

1.66 

1.8 

10 

1.086 

1.39 

9.982 

1.4 

1.015 

1.11 

1.86 

1.7 

11 

1.120 

1.45 

0.934 

1.2 

0.942 

0.826 

2.42 

1.4 

12 

1.098 

1.50 

1.13 

2.0 

0.856 

0.811 

2.90 

2.5 

standard  deviation  of  the  experimental  point  from  the  calculated  r  vs  t 
curve.  These  values  illustrate  the  scatter  of  the  data. 
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TABLE  2 

INSTANTANEOUS  LINEAR  BURNING  RATES, 

"T 

cm/sec  X  10^ 


Time  , 
seconds 

Block 

Flow 

10  .  22  Runs _ 1 

2 

3 

4 

5 

n 

t 

8 

9 

11 

12 

Average 

10 

13 

.  1 

11 

4 

11 

9 

13 

.6 

15.5 

15.8 

15.6 

15 

8 

16 

7 

18.7 

21 

5 

— 

15 

11 

1 

10 

0 

10 

1 

11 

.2 

12.5 

12.7 

12.6 

12 

9 

13 

5 

14.8 

16 

8 

12.6 

30 

8 

4 

8 

0 

7 

9 

8 

.3 

8,8 

8.9 

9.1 

9 

3 

9 

7 

10,0 

11 

3 

9.1 

45 

7 

.3 

7 

2 

6 

9 

7 

.1 

7.2 

7.4 

7.5 

7 

,8 

8 

.0 

8.1 

9 

1 

7.6 

60 

6 

.6 

6 

6 

6 

3 

6 

.3 

6.4 

6.5 

6.7 

6 

9 

7 

1 

7.0 

7 

8 

6.7 

90 

5 

8 

5 

9 

5 

.6 

5 

5 

5.3 

5.4 

5.7 

5 

9 

6 

0 

6 

7 

5.8 

120 

5 

.3 

5 

5 

5 

2 

5 

0 

4.8 

4.8 

5.1 

5 

3 

5 

4 

- 

5 

5 

5.2 

150 

5 

0 

5 

il 

4 

9 

4 

.6 

4.4 

4.4 

4.7 

4 

9 

4 

9 

- 

5 

0 

4.8 

Flow 

36,000,  21 

Runs 

10 

13 

3 

11 

3 

10 

9 

12 

3 

14.1 

14.4 

14.0 

12 

5 

L2 

9 

14 

3 

— 

15 

11 

0 

9 

6 

9 

3 

10 

2 

11.5 

11.7 

11.5 

10 

6 

ro 

9 

12 

1 

10.9 

?o 

8, 

.2 

7, 

5 

7. 

.2 

7. 

.6 

8.2 

8.4 

8 

2 

8 

.4 

8.5 

9 

.2 

8.2 

45 

6, 

.9 

6, 

5 

6, 

.2 

6, 

.5 

6.9 

7.1 

7 

2 

a 

,2 

7.4 

8, 

.0 

7.0 

60 

6. 

,2 

6. 

0 

5, 

,7 

5. 

^8 

6.1 

6.2 

6.3 

6 

5 

6 

,6 

6.7 

7. 

,2 

6.3 

90 

5. 

.3 

5. 

3 

5, 

.0 

5. 

,  1 

3.2 

5.2 

5.4 

5 

8 

5 

.8 

5.9 

6. 

.3 

5-5 

120 

4. 

9 

4. 

8 

4, 

,6 

4, 

,6 

4 . 6 

4.7 

4.9 

5 

3 

5.4 

5 

8l 

1 

5.0 

180 

4. 

3 

4. 

3 

4, 

.1 

1 

1  4. 

0 

3.9 

i 

4.2 

4 

7 

4.8 

5. 

.1 

4.4 

Flow 

6000, 

4  Runs 

60-180-" 

3. 

8 

3. 

6 

3 

3 

5 

B 

B 

m 

3 

8 

B 

D 

00 

•Constant  with  time 
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Figure  3«  Instmtaneous  Regression  Rates  at  Constant  Time 
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time 

f  (av) 

60 

4.2 

80 

4.0 

100 

3.9 

180 

3.9 

This  variation  is  due  to  the  changes  in  rates  prior  to  60  seconds  and  is 
slight,  but  does  illustrate  how  misconceptions  may  evolve  from  averages. 

There  are  two  other  ways  of  presenting  the  data:  as  linear  rates  at  con¬ 
stant  radius  (Table  3  );  and  as  total  mass  burning  rates  (Table  4)* 

The  numbers  for  Table  4  are  obtained  through  the  use  of  £q.  1  ,  2 

and  3  ,  and  are  expressed  per  cm  of  axial  length; 


(3) 

M  =  2p  TT  rr  (for  instantaneous  f)^ 
p  =  density  of  polymer 


M  =  jt)  ff  -  r^)/Lt  =  p  ir  +  rj)(r2  -  rj)/At 

where 

't  =  1  cm 


It  can  be  seen  from  “^able  4  that  the  changes  in  M  are  very  small  after 
60  seconds  . 


lim  (pg  -  rj)/A  t  =  r  and  (r^  +  r^)  =  =  2r  for  small  Ar's. 

A  t  -♦  0 
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INSTANTANEOUS  LINEAR  BURNING  RATES  AT  CONSTANT  RADIUS, 
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An  attempt  to  run  the  combustor  at  higher  pressures  (about  30  psig)  was 
not  successful  because  of  vibration  and  chatter  of  the  apparatus.  These 
runs  were  discontinued  because  of  safety  factors  and  of  considerable  doubt 
as  to  the  significance  of  the  data  obtainable  under  these  conditions.  Mod¬ 
ifications  of  the  apparatus  would  be  necessary  before  pressurized  runs  could 
be  made. 


TEMPERATURE  PROFILE  OF  RIEL  BL0(XS 

It  was  found  that  the  temperature  profile  of  the  fuel  blocks  behind  the 
regressing  surface  followed  Eq.  4. 

T  =  exp  (-bX)  (4) 

where 

T  3  temperature  above  initial  block  temperature  at  X  centi¬ 
meters  from  the  surface 

T^  -  surface  temperature  above  initial  block  temperature 

The  derivation  of  this  equation  is  discussed  in  a  following  section. 
Table  3  summarizes  the  results  of  thermocouple  measuremenisand  Fig.  3 
indicates  a  typical  fit  of  data  and  equation.  These  values  are  consid¬ 
erably  higher  than  those  previously  reported  (Bef.  3).  The  reason 
for  the  difference  is  that  0.6-  to  1-mil  thermocouples  were  used  in  the 
later  experiments,  whereas  previously  10-mil  wire  was  used.  Heat  losses 
along  wires  as  small  as  1  mil  are  believed  to  be  negligible. 
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C 


IAfil£  3 

TEMPE&ATUBE  PR0F1I£S  OF  FUEL  BLOCKS 
T  =  exp  (-bX) 


''s- 

C® 

^>-1 

cm 

Material 

36,000 

365 

9.03 

FMMA 

330 

8.88 

FMMA 

310 

11.39 

FMMA 

325 

11.00 

FMMA 

315 

12.20 

FMMA 

360 

11.03 

FMMA  +  l/4  C 

360 

9.30 

FMMA  +  l/4  f  C 

6000* 

350 

6.53 

FMMA 

370 

6.90 

FMMA 

Average  345  ^ 

’  20 

\ 

^he  temperature  profiles  were  linear  for  the  first  40  and  33  seconds, 
respectively,  before  becoming  exponential. 


c 


H-5179 


17 


A  DIVISION  OW  NOATH  AMCAICAN  AVIATION.  INC 


DISTANCE  FROM  SURFACE  X.  CENTIMETER 


Figure  Tenroeratore  of  Fuel  Behind  the  Surface 


A  OlvisioN  or  NOirrM  ambrican  aviation,  ino 


Two  important  observations  should  be  noted.  First,  the  surface  tempera* 
tures  are  the  same  (within  experimental  deviation)  for  opaque  (carbonized) 
and  clear  blocks.  Thus,  radiation  does  not  appear  to  influence  these 
temperature  readings.  Secondly,  the  temperature  profiles  obtained  with 
the  lower  flowrate  (6000)  were  almost  linear  for  about  the  first  half  of 
the  run  time.  After  this  initial  period,  they  followed  the  exponential 
form.  The  initial  temperature  rise  registered  by  the  thermocouple  was 
more  rapid  for  the  6000  flowrate  runs.  This  was  probably  due  to  smaller 
burning  rates,  and,  thus,  less  heat  was  used  up  in  pyrolysis. 

It  is  apparent  that  small  changes  in  surface  temperatures  cannot  be  de¬ 
tected  by  this  technique  because  of  the  relatively  large  scatter.  The 
deviations  are  probably  due  to  uncertainties  in  the  exact  position  of 
the  thermocouple  with  respect  to  the  surface. 

Attempts  were  made  to  obtain  optical  p3rrometer  measurements  of  the  flame 
through  the  FMMA  walls.  It  was  considered  that,  although  the  absolute 
value  of  the  temperatures  would  be  inaccurate,  trends  with  time  and  axial 
distance  may  be  meaningful.  The  values  ranged  from  1080  to  1340  K,  but 
too  much  scatter  was  obtained  for  interpretation. 


11-5179 


19 


FORM  ROt  •  FLATC  RIV.  !••• 


A  DIVISION  or  HOPTTM  AMKRICAN  AVIATION*  INO 


DISCUSSION 


Discussion  of  the  experimental  results  has  been  divided  into  two  sections. 
The  first  is  concerned  with  the  surface  temperature  measurements  and  the 
isothermal  pyrolysis  rate  data  that  are  available.  The  second  is  con¬ 
cerned  with  the  burning  rate  data  obtained  in  the  combustor^  and  compari¬ 
son  of  these  data  with  the  presently  accepted  theory.  These  sections  also 
contain  tabulations  and  discussions  of  data  considered  only  in  a  qualita¬ 
tive  manner  in  previous  progress  reports. 


PYROLYSIS  VS  OHDATIVE  DEGRADATION 

Temperature  Measurements 

Derivation  of  the  Heat  Conduction  Equation.  Equation  4  was  obtained  from 
Eq.  5,  which  is  the  differential  equation  for  a  moving  boundary  cylindrical 
surface  heat  conduction  problem  for  a  semi-infinite  solid: 


dT/dt  =  kd^T/d  +  [ (k/a)  +  r]  dl/dX  =  0 


(5) 


or 

d*^T/dX^  =  -b  dT/dX 


where 

k  ==  thermal  diffusivity  of  the  fuel,  sq  cm/sec 

X  =  linear  disteince  of  the  thermocouple  from  the  surface 

a  =  distance  from  center  of  hole  to  thermocouple,  a 
constant  (R  +  X) 
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T  =  temperature  at  X,  referenced  to  initial  temperature  of  block 

=  extrapolated  surface  temperature,  referenced  to  initial 
temperature  of  block 

b  =  (l/a  +  r/k) 

r  -  average  regression  rate 

To  arrive  at  Eq.  4,  it  vas  necessary  to  assume  that  temperatures  and  re¬ 
gression  rates  were  independent  of  time  (steady  state  and  constant 
coefficients).  The  good  fit  of  the  data  to  the  simple  exponential  form 
indicates  that  these  assumptions  were  justified  (Fig.  5)- 

The  average  value  for  T  is  345  -  20  (average  derivation)  above  initial 
temperature  (25  C),  370  C.  A  difference  of  10  centimeters  in  X  would 
change  the  extrapolated  T^  value  about  30  degrees;  therefore,  the  uncer¬ 
tainty  in  X  can  easily  account  for  the  deviations  observed. 


Correction  of  the  Extrapolated  Tg.  Because  of  the  endothermic ity  of  the 

polymer  degradation,  the  extrapolated  surface  temperature,  as  obtained 

from  Eq.  4  is  low.  The  heat  absorbed  during  reaction  gives  a  more  rapid 

drop  in  temperature  with  distance  than  would  be  predicted  from  simple  heat 

conduction  theory.  The  thickness  of  the  reacting  layer  is  X^,  at  which 

distance  the  temperature  has  dropped  to  T  .  The  temperatures  at  all 

r 

values  of  X  greater  than  X^  will  then  follow  Eq.  4.  If  it  is  assumed 
that  the  pyrolysis  rate  at  T^  is  about  100  times  slower  than  the  rate  at 
T^,  the  actual  surface  temperature,  then  T^  can  be  calculated  for  each  X^ 
using  Eq.  6  obtained  from  chemical  kinetics: 

/K^  =  100  =  exp  C-E/R(l/T^  -  1/T^)]  (6) 

a  r 
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where 

K  -  rate  constaDts 

E  -  activation  energy  -  33,000  (Ref.  6) 
and  are  in  degrees  K 
R  ~  the  gas  constant  -  1.99 


The  T  's  are  obtained  from  Eq.  4  for  each  X  .  The  results  of  these 
*  r 

calculations  are  shown,  in  Table  6,  using  approximately  the  average  ex¬ 
perimental  values  in  Eq.  4  and  =  36O  C. 

TABLE  6 

COHHECTION  OF  SURFACE  TEMPERATURE 
(360  C) 


Tr. 

c 

T 

^a 

C’ 

Correction, 

C 

35 

260 

360 

0 

3.0 

275 

15 

2.0 

300 

50 

1.0 

330 

450 

90 

0.5 

345 

475 

115 

0.1 

355 

490 

130 
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The  largest  value  of  0.035  centimeter  for  represents  the  conduction 

profile,  since  T  =  T  =  360  C,  which  would  be  obtained  only  if  a  negli- 
a  8 

gible  fraction  of  the  heat  was  absorbed  for  pyrolysis,  i.e.,  heat  conduction  > 
rate  of  heat  absorbed.  Thus,  it  is  the  maximum  value  possible  for  an  endo¬ 
thermic  pyrolysis.  The  smallest  value  tabulated,  10  ^  centimeters,  is  the 
minimum  value  to  be  examined  because  the  correction  will  not  change  signi¬ 
ficantly  for  smaller  distances. 


Comparison  of  IsotheriLal  Pyrolysis  Rates  to 
Regression  Rates 

The  next  step  is  to  compare  the  rate  of  pyrolysis  at  the  corrected 

temperature,  T  ,  as  predicted  from  isothermal  pyrolysis  rate  data,  to  the 
a 

observed  regression  rates  and  temperatures.  Presented  below  is  a  tabulation 
of  the  PMM/l  pyrolysis  rate  constants  measured  in  vacuum  by  Madorsky 
(Ref.  6),  which  followed  Eq.  7. 

K  =  3.2  X  10^  exp  (-16,500/t)  (min“^)  (?) 


Temperature, 

K  X  10 
.  -1 

C 

min 

240 

1.9 

250 

3.7 

260 

6.9 

270 

11.0 

These  values  were  obtained  by  extrapolating  the  rate  data  to  0^  decom¬ 
posed,  and  represent  the  upper  limit  of  the  data. 
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On  re-examination  of  the  pyrolysis  rate  data  obtained  at  Rocketdyne, 
it  was  found  that  some  values  reported  (Fig.  11^  Ref.  2)  were  in  error. 

The  corrected  values  are  presented  and  discussed  in  Appendix  A,  There 
is  general  agreement  between  the  data  obtained  by  Madorsky  (Ref,  6)  and 
that  obtained  in  this  laboratory. 

It  is  now  necessary  to  take  the  corrected  T^,  or  T^,  and  the  rate  para¬ 
meters  from  Eq.  7  and  calculate  the  theoretical  mass  burning  rates.  First, 
a  rate  equation  to  be  used  in  a  system  with  an  exponential  temperature 
gradient  must  be  developed.  The  development  is  presented  below. 

For  a  first  order  reaction,  such  as  PMMA.  pyrolysis  (Ref.  6),  the  mass 
pyrolysis  rate  is  proportional  to  the  mass  heated  and  the  rate  constant, 

Eq.  8. 


M  =  Mf  exp  (-E/RT) 


(8) 


where 


ft  =  mass  pyrolysis  rate 

M  =  mass  heated  to  temperature,  T,  which  is  in  degrees  K 
f  =  pre-exponential  factor 
E  =  activation  energy 
R  =  gas  constant 

If  the  mass  pyrolysis  rate  per  unit  area,  m,  is  considered,  M  can  be  re¬ 
placed  by  p  V/A  (density  times  volume  divided  by  surface  area),  which  re¬ 
sults  in  pX  (where  X  is  the  distance  from  the  surface).  Thus,  Eq.  9  can 
be  written; 

m  =  p  Xf  exp  (-E/RT)  (q\ 
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Since  the  temperature  is  a  function  of  X,  the  total  m  is  the  sum  of  the 
rates  over  the  entire  range  of  temperatures;  thus,  Eq.  10  results: 

8 

m  =  p  f  exp  (-E/RT)dX  (lO) 

where 


s  -  thickness  of  the  polymer. 

The  temperature  dependence  on  X  as  expressed  by  Eq.  11  can  be  substituted 
for  T: 


t'  =  t'  exp  (-cX)  (ll) 

a 

Temperatures  (T  )  must  be  in  K,  and  can  be  calculated  from  the  values  for 

X  ,  T  ,  and  T  in  Table  6  The  difference  between  Eq.  4  €tnd  11  is  the  tem^ 
r  r  a 

perature  scale  and  the  exponential  coefficients  in  the  latter  p3rrolysis  is  considered. 
Equation  11  is  a  simple  exponential  which  approximates  satisfactorily  Eq.  4 
for  small  values  of  X,  which  is  the  region  of  importance  for  the  pyrolysis. 

If  a  new  variable,  y,  is  defined  by  Eq.  12:. 


y  =  E/RT  =  (e/RT  )  exp  (cX) 
a 

(12) 

and 

dy  =  (cE/RT  )  exp  (cX)dX 
a 

or 

dX  =  dy/cy 

then  Eq.  10  through  12  can  be  combined  to  give  Eq.  13: 


m  ■=  p  f/c  J  exp  (-y)/y  dy  (l3) 

e/rt^ 
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The  upper  limit  of  ®  can  be  used  in  place  of  a,  since  the  value  of  the 

integral  is  not  changed  significantly  at  values  of  X  larger  than  X  be- 

r 

cause  of  the  rapid  decrease  in  temperature  (the  error  introduced  by  using 
®  is  about  1  percent) .  The  value  of  this  integral  has  been  evaluated  by 
a  series  solution.  The  results  of  the  calculated  and  observed  burning 
rates  are  presented  in  Table  7  . 


TABLE  7 

CALCULATED  MASS  BURNING  RATES 


Tr, 

K 

K 

Xj.  X  100, 

centime ters 

c 

00 

J  «  10*’ 

e/rt^ 

P  f/c  X  10® 

m  X  10  , 
gm/sec/sq  cm 

533 

633 

3.5 

4.9 

1.8 

7  8 

1.4 

548 

648 

3.0 

5.6 

3,3 

6.8 

2.2 

573 

683 

2.0 

8.7 

12.8 

4.4 

5.7 

603 

723 

1.0 

18.1 

52 

2.1 

10.9 

618 

748 

0.5 

38.2 

115 

1.00 

11.5 

628 

763 

0.1 

^95 

180 

j  0.20 

3.6 

Experimental*  m  =  100  x  10  ^  gm/se</sq  cm 


^This  value  represents  an  average  number  over  about  10  runs  of  60  second 
duration,  24-inch  axial  length,  and  36,000  flov. 

It  can  be  seen  that  the  predicted  are  too  lov  by  a  factor  of  about  10, 
even  using  the  highest  rates  calculated  from  the  highest  isothermal  rates. 
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Some  explanations  that  can  be  advanced  at  this  time  for  the  low  calculated 
values  are: 

1.  The  polymer  regression  rate  is  not  governed  solely  by  rolysis; 
oxidative  degradation  contributes.  This  may  be  caused  by  O2, 

CO^,  CO,  H^O  or  various  radical  intermediates  initiating  chain 
degradations , 

2.  Some  fuel  is  blown  off  the  surface  as  the  polymer  softens.  This 
appears  unlikely  for  PMMA  but  may  be  significant  in  the  higher 
rates  observed  for  polystyrene  (PS)  since  it  softens  at  lower 
temperatures . 

The  temperature  measurements  are  significantly  low  due  to  heat 

losses  An  approximate  calculation  shows  that  if  the  temperature 

correction  is  in  error  by  80  K  and  the  X  is  about  0.02  cm,  the 

r  ^ 

calculated  m  is  still  too  low  by  a  factor  of  2  (m  -  49  x  10  ). 

4.  The  extrapolation  of  the  isothermal  rate  data  may  be  in  error  and 
the  rates  are  maiich  faster  than  Eq,  7  predicts. 

It  is  believed  that  the  above  calculation  reveals  a  serious  discrepancy 
and  that  illucidation  of  these  differences  may  lead  to  a  much  better  under¬ 
standing  of  hybrid  combustion.  The  ability  of  oxygen  to  increase  the  rate 
of  depolymerization  of  PMMA  by  as  much  as  a  factor  of  10  is  shown  by  the 
resuits  tabulated  in  the  Appendix  Table 
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COMBUSTION  CONSIDERATIONS 

The  combustion  rate  data  which  were  obtained  during  the  course  of  the 
program  are  discussed  in  this  section*  Most  of  the  experiments  were 
conducted  under  oxidizer-rich  conditions*  For  FMMA  to  burn  to  CO  or 
COg  and  water,  it  is  necessary  to  have  mixture  ratios  about 

1*1  and  1.9,  respectively;  as  can  be  seen  from  Table  9  and  10,  the 
ratio  employed  experimentally  is  usually  greater* 


Instantaneous  Burning  Rate  Behavior 

Table  2,  3,  and  8  summarize  the  instantaneous  burning  rate  behavior* 
The  calculations  of  slopes  (r  vs  t)  at  the  extreme  times  (lO  and  150 
or  180  seconds)  are  not  too  reliable  because  of  a  lack  of  data  beyond 
these  points* 

General  trends  can  be  detected.  A  slight  minimum  occurs  for  the  two 
higher  flowrat€.*s  and  appears  to  shift  to  longer  axial  distances  at 
larger  radii  for  the  55,000  flow.  Table  This  minimum  is  probably 
still  caused  by  an  entrance  effect  since  it  is  much  less  pronounced 
with  the  smoother  entrance  section.  Block  1  is  not  tabulated  in  most 
cases  since  its  behavior  is  anomolous  because  of  the  pronounced 
entrance  effect. 

The  linear  regression  rates  decrease  with  time,  as  predicted,  because  of 
the  larger  port  area*  The  rate  appears  to  show  an  increase  with  down¬ 
stream  distance*  At  the  smaller  radii  the  ratio  of  the  instantaneous 
regression  rates  for  Re  55,000/36,000  increases  with  axial  distance  as 
presented  in  the  previous  quarterly  report  (Ref*  5);  but  this  does  not 
hold  at  larger  radii  (Table  8).  The  fact  that  the  ratio  of  rates  remains 
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constant  or  drops  slightly  with  axial  distance »  as  is  the  case  for 
larger  radii,  would  be  predicted  from  flow  velocity  considerations, 
because  the  ratio  of  velocities  also  drops  slightly.  However,  the  rise 
observed  at  smaller  radii  cannot  be  explained  from  flow  considerations 
alone.  Similarly,  the  plateaus  and  drops  evidenced  at  the  shorter 
times  or  smaller  radii  (Fig.  3  and  4),  appear  real  and  would  require 
reversal  of  trends  in  terms  of  flow  considerations.  Other  factors 
must  be  important. 


Another  interesting  result  is  that  the  total  mass  burning  rate,  ft,  is 
relatively  constant  with  time  after  about  60  seconds.  This  is  of  con¬ 
siderable  importance  in  hybrid  grain  design.  In  the  following  section 
the  accepted  regression  rate  equation,  which  has  been  the  basis  of  grain 
designs  is  compared  to  the  data  obtained  in  this  program. 


Theory  vs  Experiment 


The  accepted  equation  for  the  average  burning  rate  as  derived  by  UTC 
(Ref.  7)  and  used  by  Lockheed  (Ref.  8)  for  grain  design  problems 
is : 


r  = 


(A)"- 


+  b' 


(14) 
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where 


a  and  b*=  complex  parameters  involved  in  heat  transfer 
(not  constants) 

w^  =  total  flowrate  of  propellants 
L  =  grain  length 
A  =  port  area 


The  b* term  is  the  contribution  of  radiation*  It  is  usually  neglected 
because  it  is  small  and  difficult  to  evaluate.  As  indicated  by  the 
present  study,  radiation  effects  on  rate  do  not  appear  important 
(Table  ll)*  Thus,  the  b*term  has  been  neglected  in  the  present  cal¬ 
culations  • 


Table  9  &nd  10  summarize  the  data  (average  f)  for  runs  employing  blocks 
of  24  inches  total  length,  at  various  times  and  flowrates.  Equation  14 
was  used  to  calculate  the  parameter  the  values  listed  indicate  the 
variation  of  ^  with  time  and  oxygen  flowrate.  Involved  in  this  parameter 
are  ratios  of  velocities  and  Stanton  numbers,  enthalpies  of  gas  and  solid 
phases,  viscosities,  etc.  In  view  of  the  complexity  of  a,  these  values 
appear  remarkably  constant  under  the  range  of  conditions  considered  in  this 
program.  Two  comments  are  warranted  at  this  time:  (l)  the  function  ex¬ 
pressed  by  Eq.  14  (a  =  approximately  constant)  appears  satisfactory  for  most 
applications,  although  the  influence  of  possible  deviations  of  a  with 
time,  flowrates,  and  especially  with  different  propellant  systems 
should  be  checked  for  their  effect  on  grain  design,  etc. 
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TABLE  9 


TIME  DEPENDENCE  OF  HEAT  IBANSFER  PARAMETER 


Time , 
seconds 

t  X  100 , 
av  ^  * 

gm/sec 

Total  Propellant 
Flow,  w^, 

gm/sec 

Port  Area ,  A, 
sq  cm 

a  X  1000 

36,000  = 

14.4  gm/sec  of 

15 

10.9 

21.7 

6.60 

1.84 

30 

8.2 

20.4 

7.95 

1.67 

45 

7.0 

19.8 

9.07 

1.62 

60 

6.3 

19.6 

10.3 

1.67 

90 

5.5 

19.3 

12.5 

1.71 

120 

5.0 

19.2 

14.7 

1.76 

180 

4.4 

19.1 

19.0 

1.93 

55,000  =  22.0  gm/sec  of  0„ 


15 

12.6 

30.7 

7.16 

30 

9.1 

28.8 

8.70 

45 

7.6 

28.0 

10.1 

60 

6.7 

27.7 

11.4 

90 

5.8 

27.3 

13.6 

120 

5.2 

27.1 

15.9 

150 

4.8 

27.0 

18.3 
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^^Initial  ID  of  l/2  inch. 
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(2)  It  is  believed  that  experimental  values  of  the  parameters  making  up 
a  are  not  sufficiently  well  known,  under  the  conditions  of  hybrid  combus¬ 
tion,  to  state  that  theory  and  experiment  are  in  agreement  (or  disagree¬ 
ment).  The  apparent  agreement  may  arise  from  lack  of  sensitivity  of  the 
rate  to  the  parameters  involved  or  from  compensating  deviations* 

Thus,  it  appears  that  too  little  is  known  at  this  time  for  a  rigorous 
combustion  model  to  be  developed;  the  need  for  more  basic  information 
concerning  both  chemical  and  physical  processes  is  still  indicated* 


Regression  Rates  of  Different  Fuels 


Although  the  regression  rate  results,  employing  different  fuels,  have  been 
discussed  in  previous  reports  (Ref»  3  and  4),  Table  11  presents  a  more 
complete  summary  of  quantitative  values. 

Certain  points  can  be  mentioned  briefly.  Preheating  or  precooling  of  the 
fuel  blocks  had  little  effect  on  the  rates*  The  addition  of  carbon  to 
PMMA  appeared  to  yield  a  small  increase  in  rate.  However,  the  10^  change 
observed  is  probably  due  to  the  difference  in  degree  of  polymerization  of 
the  fuel  in  the  presence  of  carbon,  rather  than  to  an  increase  in  heat 
transfer  from  radiation,  because  PMMA  blocks  containing  more  carbon 
than  1/2^  did  not  cure  completely.  The  small  increase  in  rate  from  PMMA 
to  30^  A1  in  PMMA  is  surprising  since  several  factors  would  appear  to 
promote  a  higher  rate,  for  example:  The  higher  flatme  temperature  of  an 
aluminized  system,  the  smaller  amount  of  polymer  to  pyrolize  per  centimeter 
of  distance  (the  aluminum  would  be  assumed  to  melt  but  not  vaporize 
between  the  wall  and  the  flame),  and  the  higher  absorbtivity  of  the 
opaque  wall.  The  60^  A1  and  30/i  TiHj^  fuels  do  show  a  large  increase  in 
rate*  Again  the  contribution  of  the  additive  to  the  degree  of  polymerization 
of  the  PMMA  is  unknown,  which  may  account,  at  least  in  part,  for  differences 
in  rates* 
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TABI£  11 

COMPARISON  OF  SURFACE  REGRESSION  RATES 
OF  FUELS 


Fuel 

r  X  1000^^ 
cm/sec 

Initial 

Fuel 

Temperature , 
C 

Og  Flow 
55,000 

6000 

PMMA 

1 

4.2 

Ambient 

PMMA.  +  30  A1 

1 

4.7 

Ambient 

PMMA  +  60  A1 

1 

7.8 

Ambient 

PMMA  +  30  ^  TiH. 

14.8 

5.0 

Ambient 

PE 

11.4 

4.8 

Ambient 

PS 

25.0 

9.8 

Ambient 

PMMA 

11.3 

- 

-20 

PMMA 

12.1 

- 

100 

A* 

cm/sec 

Fuel 

. 

PMMA 

PMMA 

Run  Time, 

PMMA 

+  1/4  56  c 

+  1/2  56  c 

seconds 

Og  Flow 

10.1 

11.2 

30 

36,000 

8.8 

9.3 

- 

53 

8.5 

CM 

- 

60 

7.6 

8.0 

8,3 

90 

9.2 

- 

10.1 

45 

^<Average  rates  over  two  60-second  runs  of  24-inch  axial  length;  older 
24 -inch  entrance  duct  used. 

■^^Average  rates,  fuels  all  at  ambient  temperature;  newer  1-inch  entrance 
section  used. 
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The  behavior  of  polyethylene  (P£)  was  similar  to  FHKA.  except  for  a  sig¬ 
nificant  amount  of  melting  and  liquid  ejection.  The  burning  of  polysty¬ 
rene  (PS)  was  considerably  different  from  that  of  the  other  fuels.  As 
discussed  in  previous  reports  (Ref.  3  end  4  )  the  rate  became  more 
rapid  immediately  following  the  entrance,  and  maintained  a  very  high  rate 
of  regression  which  started  dropping  after  about  24  inches  (35 9 000  oxygen 
flowrate) . 

A  comparison  of  the  pyrolysis  rates  of  FMMA  and  PS  (Ref.  6  )  indicates 

that  the  PS  rate  is  slower  in  the  temperature  range  studied,  contrary  to 
regression  rate  results.  An  important  contributor  to  the  observed  higher 
rate  may  be  the  lower  softening  and  melting  points  of  PS,  The  roughness 
of  the  surface  after  burning  may  be  indicative  of  the  importance  of  this 
factor.  In  addition,  the  flame  temperature  of  the  PS  oxygen  system  is 
higher,  and,  also,  if  the  monomer  pyrolizes  at  or  near  the  surface,  addi¬ 
tional  heat  will  be  liberated  to  increase  the  voltalization  rate  because 
styrene  has  a  positive  heat  of  formation. 

These  results  point  to  a  need  for  a  better  understanding  of  relative  py¬ 
rolysis  and  oxidative  pyrolysis  rates,  and  the  influence  of  additives 
(metals,  etc.,  on  these  rates).  In  addition,  the  pyrolysis  rates  of  the 
monomers  (e.g.,  styrene  and  MMA)  must  be  known  better  to  evaluate  their 
contribution  to  surface  regression  rates. 
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SUMMARY  AND  RECOMMENDATIONS 


The  burning  of  a  solid  fuel  and  a  gaseous  oxidizer  has  been  studied  using 
primarily  plexiglass  (pmma)  as  the  fuel,  and  oxygen,  at  various  flowrates, 
as  t^e  oxidizer.  The  burning  rates  were  determined  at  different  values 
of  axial  distance  and  time  to  define  experimentally  the  rates  as  functions 
of  these  parameters.  Additional  measurements  were  made  of  the  temperature 
of  the  fuel  behind  the  regressing  surface  to  determine  the  surface 
temperature , 

The  extrapolated  surface  temperatures  appear  to  be  too  low  to  account 
for  the  observed  regression  rates  by  only  a  pyrolysis  mechanism  if  the 
isothermal  pyrolysis  rates  of  PMMA  are  extrapolated  to  the  surface 
ten^erature . .  One  or  more  of  the  following  factors  may  contribute  to  the 
observed  difference:  (l)  oxygen,  radical  intermediates,  or  products  may 
accelerate  pyrolysis,  (2)  the  temperature  measurements  may  be  in  error, 
and  (3)  the  extrapolated  pyrolysis  rates  may  be  in  error.  The  reason  for 
this  difference  must  be  clarified.  If  (l)  is  important  then  extrapola¬ 
tion  of  behavior  to  different  propellant  systems  will  be  extremely  dif¬ 
ficult  unless  the  interactions  and  relative  reaction  rates  are  known. 

Similarly,  the  high  regression  rates  of  polystyrene  (PS),  greater  than 
two  times  that  of  PMMA,  can  be  explained  only  in  part  by  the  difference 
in  flame  temperature.  The  heat  of  depolymerization  plus  vaporization  of 
the  monomer  of  PS  and  PMMA  are  about  263  and  208  cal/gm,  respectively. 

Thus,  the  heat  transfer  rate  necessary  for  the  PS  burning  rate  would  be 
almost  three  times  that  required  for  PMMA,  which  appears  unlikely.  It 
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is  possible  that  the  pyrolysis  of  styrene  monomer  could  account  for  addi¬ 
tional  heat  since  its  heat  of  formation  is  about  (+)  35  kcal/mole 
(endothermic),  but  the  pyrolysis  must  occur  sufficiently  close  to  or  on 
the  surface  (thus  rapidly)  to  yield  the  desired  effect.  Again,  the  role 
of  chemical  kinetics  must  be  explored  further  because  heat  transfer  and 
thermochemistry  apparently  do  not  account  for  the  observed  results  under 
the  conditions  of  this  program. 

In  conclusion,  the  heat  transfer  considerations  used  to  explain  regression 
rate  behavior  are  important.  There  are  indications,  however,  that  chemical 
kinetics  may  also  contribute  significantly  to  burning  rates  for  hybrids, 
and  that  these  contributions  appear  more  pronounced  in  the  comparisons 
between  different  propellant  combinations. 
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APPENDIX  A 


This  appendix  presents  tabulations  of  the  rate  constants  determined  for 
PMMA  pyrolysis. 

On  thorough  re-examination  of  the  data  reported  previously  (Ref.  2  Fig.  11  ) 
it  was  found  that  values  of  two  points  were  in  error  in  the  original  cursory 
calculations;  the  corrected  results  are  presented  in  Table  12  .  It  is  be¬ 
lieved  that  the  results  agree  fairly  well  (within  a  factor  of  two)  with 
those  reported  by  NBS  (Ref.  6  ) (using  their  average  rates  rather  than 
those  obtained  by  extrapolating  to  0  ^  decomposed),  considering  the  large 
difference  in  size  of  the  polymer  saii^)les  (a  factor  of  about  250). 

The  results  in  Table  12  are  the  average  rates  over  their  entire  range  of 
pyrolysis,  which  appear  to  be  the  best  representation  of  the  data.  No  con¬ 
sistent  trends  with  percent  volatilized  were  noted  in  these  data. 

Figure  6  and  7  illustrate  Arrhenious  plots  of  the  data  as  first-  and 
zero-order  pyrolyses  reactions,  respectively.  It  appear within  the 
limitations  of  the  accuracy  of  the  experiments,  that  the  data  can  be  re¬ 
presented  almost  equally  well  by  either  order. 

It  is  believed  that  the  rates  reported  at  the  higher  temperatures  (Ref.  2  ) 
gave  nonlinear  Arrhenious  plots  because  of  the  limitations  of  the  experi¬ 
mental  techniques  used,  rather  than  because  of  changes  in  chemical  mechanism. 
For  this  reason,  they  have  not  been  used  in  the  present  report.  This  does 
not  imply  that  changes  in  mechanism  do  not  take  place  at  higher  temperatures 
but,  rather,  that  the  higher  rates  measured  by  the  techniques  used  in  this 
program  cannot  be  interpreted. 
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TABLE  12 

FMMA  mOlAfSIS  RATES 


Gas 

Temperature , 

C 

10  Vt 

Initial 

Weight, 

grams 

Average 

Bate  (kj, 
mg/min 

k,  X  10^ 
1>  .-1 
min  ^ 

^2 

256 

1.89 

1.33 

1.5 

1.3 

^2 

269 

1.84 

1.57 

2.0 

1.4 

^2 

318 

1.69 

1.42 

5.4 

4.1 

2^0^ 

230 

1.99 

1.37 

2.0 

1.6 

2i>\ 

242 

1.94 

1.48 

3.4 

2.5 

2  %  Og 

263 

1.87 

1.55 

4.6 

3.5 

2  $  Og 

281 

1.80 

1.52 

6.0 

4.7 

2  ^  Og 

312 

1.71 

1.58 

11.3 

10.8 

35  ^  Og 

210 

2.07 

1.58 

1.1 

0.8 

35  Og 

225 

2.00 

1.67 

3.2 

2.0 

35  ^  Og 

235 

1.97 

1.44 

3.4 

2,6 

35  ^  O2 

249 

1.92 

1.22 

6.0 

5.5 

35  ^  Og 

262 

1.87 

1.82 

12.3 

8.7 

100  )t  Og 

203 

2.11 

1.27 

0.9 

0.7 

100  Jl  Og 

235 

1.97 

1.56 

5.1 

3.6 

CM 

0 

0 

0 

s-» 

247 

1.92 

1.66 

5.7 

4.1 

»-• 

0 

0 

0 

10 

263 

1.87 

1.69 

10.0 

7.0 
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Figure  6«  First -Order  Rate  Constants  for  FMM^  Pyrolysis 


Rq  mg /min 
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Figure  7*  Zero-Order  Bate  CoMtt«te  for  ^rrolyaie 
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